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way of interpolation is also reflected in the orthogonal polynomials corresponding to 
such a random matrix ensemble. Although the chiral Gaussian unitary ensemble 
as well as the Gaussian unitary ensemble are associated to the Dyson index /? = 2 
the intermediate ensembles exhibit a mixing of orthogonal polynomials and skew- 
orthogonal polynomials. We consider the Hermitian as well as the non-Hermitian 
Wilson random matrix and derive the corresponding polynomials, their recursion 
relations, Christoffel-Darboux-like formulas, Rodrigues formulas and representations 
as random matrix averages in a unifying way. With help of these results we 
derive the unquenched /c-point correlation function of the Hermitian and then non- 
Hermitian Wilson random matrix in terms of two flavour partition functions only. This 
representation is due to a Pfafhan factorization drastically simplifying the expressions 
for numerical applications. It also serves as a good starting point for studying the 
Wilson-Dirac operator in the e-regime of lattice quantum chromodynamics. 
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1. Introduction 



In the microscopic limit chiral random matrix theory (xRMT) can be directly mapped 
to the e-regime of quantum chromodynamics (QCD) and is successfully applied to it 
since the 90 's [HE]- Both theories share the same universality class which is the reason 
for the existence of this equivalence. xRMT was also extended to a non-zero chemical 
potential by adding a scalar proportional to 70 [S IH El [6], [7] . In the last decade a second 
appoach was pursued. A second chiral random matrix was introduced yielding the chiral 
analogue of the Ginibre ensembles [Bl El [IDl IHl 1121 [13, IIH US] • A quantitative analysis 
of the sign problem in Monte-Carlo simulations was quite elusive until it was solved in 
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xRMT [HI [H Uni |T7]. The hope is now to extend these new insights to QCD at non-zero 
lattice spacing. 

Recently, random matrix theories for lattice QCD became the focus of interest. 
The idea is to derive analytical results of lattice artefacts in the data. One important 
realization of lattice QCD is by means of staggered fermions. In Refs. IT9] . a xRMT 
was considered which is equivalent to the e-regime of these fermions. Unfortunately, this 
model is highly involved due to the high number of low energy constants and, hence, of 
the coupling constants in the random matrix model. 

The Wilson Dirac operator is another realization of lattice QCD. It proved that the 
corresponding random matrix model [201 ED [22l EHl |T9l El] is much better accessible 
for analytical calculations than the one of the staggered fermions. The Wilson term 
which is given by a Laplace operator [25| 126] explicitly breaks chiral symmetry and is 
Hermitian. Thus the main idea was to add on the diagonal of xRMT two Hermitian 
matrices to simulate the same effect [20] and it proved to be in the same universality 
class as the Wilson Dirac operator in the e-regime [271 EH ESI ED] . Actually one can 
consider a Hermitian version [201 ED E21 E31 ED E2] of this random matrix ensemble 
which is numerically cheaper in lattice simulations. However only the non-Hermitian 
version [191 EH E2] is directly related to the chiral symmetry breaking by a finite lattice 
spacing. The Hermitian version can also be considered as an interpolation between a 
chiral Gaussian unitary ensemble (xGUE) and a Gaussian unitary ensemble (CUE). 
The coupling constant is then the lattice spacing. 

Quite recently this new kind of random matrix model has given new insights on the 
signs of the low energy constants in the chiral Lagrangian of the Wilson Dirac operator 
[211 133] El]. These signs are controversial since they are crucial to decide if an Aoki 
phase [35] exists or not. Such a phase is a pure lattice artefact and has no analogue 
in continuum QCD. Therefore a large analytical [23 ESI E21 El] as well as numerical 
[371 EH ESI ED E2] effort was made to determine the low energy constants. 

Orthogonal polynomial theory [lOlilDllSlllSlllllllSlllSlIlT] was as successfully 
applied to RMT as the supersymmetry method [IHlllSlEOlEDESllSSlEllESlEe]. In 
particular the combination of both methods with the recently developed method of an 
algebraical rearrangement of the joint probability density with quotients of characteristic 
polynomials ^^7\ are quite efficient to find compact and simple analytical results of the 
spectral correlations of random matrix ensembles. In this work we address the /c-point 
correlation functions of the Hermitian as well as the non-Hermitian version and make 
use of such a combination. 

An interesting point of view of Wilson RMT appears when we study it with help 
of orthogonal polynomial theory. In Ref. [23] the authors considered the Hermitian 
Wilson RMT and found that the construction of the skew-orthogonal polynomials 
strongly depend on the index u of the random matrix which is the number of zero 
modes in the continuum limit. They only explicitly constructed these polynomials 
for u = 0,1. In this article we construct these polynomials for an arbitrary index 
and for both version of Wilson RMT in a unifying way. We also successfully look for 



Mixing in orthogonal polynomial theory 



3 



a recursion relation, Christoffel Darboux-like formulas, Rodrigues formulas and their 
explicit expression as random matrix integrals. By this study we get a complete picture 
what these orthogonal and skew-orthogonal polynomials are and how are they related 
to the orthogonal polynomials of some limits, in particular the continuum limit and the 
limit of a large lattice spacing. 

First we specify what are the conditions the orthogonal and skew-orthogonal 
polynomials have to fulfill. Thereby we recognize that they can only fulfill the conditions 
if the corresponding weight has a particular property. Luckily we are able to modify 
the weight in the joint probability densities without changing the partition functions 
and the fc-point correlation functions such that this property is fulfilled. In the second 
step we construct the polynomials with help of Pfaffians whose anti-symmetry under 
permutations of rows and columns proves quite useful. 

After we show some useful properties of the orthogonal and skew-orthogonal 
polynomials we derive a Pfaffian factorization of the fc-point correlation functions which 
is numerically advantageous because it reduces the complexity of the integrand to an 
average over two characteristic polynomials. In combination with the supersymmetry 
method [511 E21 EH ES] one may simplify the whole problem to two-fold integrals. 
Factorizations to determinants and Pfaffians were found for many random matrix 
ensembles of completely different symmetries jlOllillllSlESllllllMlESlET]. A Pfaffian 
for the eigenvalue correlations of the Hermitian Wilson random matrix ensemble was 
already shown in Ref. [23]. We also prove the existence of such a structure for the 
non- Hermitian version. Furthermore we identify the kernels of both Pfaffians with two 
flavour partition functions. The identiflcation as well as the structure carry over to 
the microscopic limit which makes them also applicable to the chiral Lagrangian of the 
Wilson-Dirac operator. 

We consider unquenched Wilson RMT, i.e. a flnite number of fermionic flavours. 
Recently, the partition function with one fermionic flavour and the corresponding 
microscopic level density was studied in Ref. [50] • In our calculations the number of 
fermions may be arbitrary. Nevertheless all eigenvalue correlations, also the one of the 
unquenched theory, can be expressed by two-flavour partition functions because of the 
Pfaffian factorization. 

Moreover such a Pfaffian determinant of the fc-point correlation functions plays 
into the hands when calculating the individual eigenvalue distributions. The authors of 
Ref. ^1] were able to express the gap probability of the eigenvalues of the Hermitian 
Wilson random matrix ensemble as a Fredholm- Pfaffian only due to this structure. 
Hence a similar simplification is highly desirable for the non-Hermitian Wilson random 
matrix ensemble. 

The outline of this article is as follows. In Sec. H] we briefly introduce the Wilson 
random matrix model and its two kinds of joint probability densities corresponding to 
the Hermitian and the non-Hermitian version. With help of the fc-point correlation 
function we propose the problem. In particular we will list the conditions the 
polynomials have to fulflU. In Sec. |3] we contruct the orthogonal and skew-orthogonal 
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polynomials. Thereby we derive a recursion relation which is helpful to proof 
a Christoffel Darboux-like formula. Moreover we show a representation of the 
polynomials and the Christoffel Darboux-like formula as random matrix averages. Such 
a representation is useful to study the microscopic limit of the random matrix ensemble 
by means of the supersymmetry method. Hence we calculate the asymptotics of the 
polynomials and the Christoffel Darboux-like formula. In Sec. H] we apply the derived 
results to the A;-point correlation functions of the Hermitian and the non-Hermitian 
version of Wilson RMT and identify the kernels of the Pfaffian with two flavor partition 
functions. Readers only interested in the /c-point correlation functions of the Wilson 
random matrix ensemble can jump to this section because it contains the main results 
which can be mostly understood without the technical details in Sec. [3] due to the 
identification of the kernels with two flavor partition functions. The conclusions are 
made in Sec. [5] and the details of the calculations are given in the appendices. 



2. Two joint probability densities for one random matrix theory 

The models we want to consider are motivated by the Wilson Dirac operator in lattice 
QCD [20]. The corresponding random matrix theory consists of the matrix 

- I 'I I (2.1) 



distributed by the Gaussian 



( 2 n + u 2 



(2.2) 



X exp 



:(trv4^ + tr5^) -ntiWW^ + fi, trA + /ii ti B 



2a2 

The Hermitian matrices A and B have the dimensions n x n and (n + z/) x (n + u) and 
explicitly break the chiral symmetry, 

75 Dw\rn=ol5 7^ - Dw\m=o "^^^^ ^5 = diag (1„, -tn+u)- (2.3) 
The matrix W is a n x [n + u) complex matrix with independent entries. The variable 
a plays the role of the lattice spacing and the Gaussian of A and B yields one low 
energy constant known as Ws [201 ED [221 [121 [211 [M]- The variables /ir/i might be also 
considered as Gaussian distributed random variables and generate two additional low 
energy constants, Wq and Wj [211 [221 [M], in chiral perturbation theory of the Wilson 
Dirac operator [271 [211 [211 [30] . Here we consider them as fixed constants to keep the 
calculation as simple as possible but the model is general enough to introduce also the 
Gaussian integrals for /ir/i at the end of the day. 

The parameter u is called the index of the Dirac operator and is the number of the 
generic real modes of -Dw Since is 75-Hermitian, i.e. Dy^ = ^^D^^^^, the matrix 

= i^D^ (2.4) 

is Hermitian. Moreover the complex eigenvalues of come in complex conjugated 
pairs only. The number of these pairs, /, varies from to n. 
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The matrix Z^s = VxV^^ can be diagonalized by a unitary matrix V eU {2n + u) 
whereas the matrix = UZiU~^ can only be quasi diagonahzed by a matrix in the 
non-compact unitary group U G U(n, n + v), i.e. = 75^7^75. The diagonal matrix 
X = diag {xi, . . . , X2n+v) consists of real eigenvalues, only. The quasi- diagonal matrix 



(2.5) 
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/ 



depends on the real diagonal matrices x^^-* 

.(2)a „(2) 



XI 'J, y'-'-' = diag {y^\ . . . , y}^^) and x'^^^ = diag {x^i', ■ ■ ■ , x)^^j^_i) with the dimensions /, 



,{2)> 



diag (x^ 

„(3) 



(1) 



(2) 



diag {x\ 



(2) 



„(3) 



I, n — I and z/ + /, respectively. Then the complex conjugated eigenvalue pairs of 
Dy/ are {z'^'^\ z*^'^'^) = (x*^^^ + ty^'^\x^'^^ — ty^"^^). The n + 1 different sectors of different 
numbers of the complex conjugated pairs are labelled by /. 

The joint probability density is one of the best quantities for analyzing the 
eigenvalue correlations of random matrices. It is also the starting point of our discussions 
in the ensuing sections. The Hermitian, D5, and the non- Hermit ian, -Dw ? Wilson 
random matrix have different joint probability densities. Though these densities have 
a completely different form, we will see that their orthogonal and skew-orthogonal 
polynomials have much in common, see Sec. [31 
The joint probability density of is [23| 

2 



P5 



{x)d[x] = c_(l - a2)-"("+'^-i/2)^-n-.^ 



exp 



X Pf 
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l<ij<2n+i/ l<i<2n+u 



{-x]-'g^^\x,)} 



where 

g''2\xi,X2) = exp 



n 



l<i<u 
l<j<2n+i' 



n 



4a2 



[X1+X2) - -{xi ~ X2) + 







nX-, 



2n+u 

n ^^j' 



(2.6) 



4a? 



-(xi + X2) 



(2.7) 



g\ '{x 



X erf 



exp 



'n{x2 - Xi) - mg-) , (ra(xi - X2) - mg-) 



n 



2a2 

We define the constants 



X ± yUlX 



a± 



2a| 

n 



2(1 ±a2)' 

(/Xr + Ail), 



2a| ^ 
n 



(2.8) 

(2.9) 

(2.10) 
(2.11) 
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IGvr 



n 



n/2 n—1 n+u—1 

(27r)'^/'n-'^'/2-n(n+.)^2n + z/)! JJ j! JJ 



(2.12) 



The notation of mgi and Xi± reflects the nature of their symmetries. The constant mgi 
acts as an effective mass and Ay^ as an effective source variable proportional to 75. They 
refer to the low-energy constants Wq and Wj which are found in the microscopic limit 
[201 ED [221 [ISl [21E1], i.e. n ^ 00, a = ^/na/V2 = a± = const., mg = mgi = const., 
A7 = = const, and £ = 2nz = const. Please notice that our notation differs from 
the one in Refs. [201 121] where the source terms proportional to 75 are denoted by z. 
To avoid confusion with the complex eigenvalues of we denote these variables by A 
in the present article. Furthermore we renamed the variables i/q and i/f to mg and A7, 
respectively, since the former notation can create a confusion with the imaginary parts 
of the complex eigenvalues of • 

The Vandermonde determinant is given by 

A2n+.(a:) = n - = (-l)(^"+^)(^"+'^-^)/^ det [xl'] . (2.13) 

l<i<j<2n+u 

The function erf(xi,X2) = erf(x2) — erf(xi) is the generalized error function. 

The Pfaffian in ps, see Eq. (12. 6p . is due to the symmetrization of the eigenvalues. 
The two point weight g2~^ is antisymmetric and is a strong interaction of two different 
eigenvalues. In the continuum limit, a — )■ 0, generates a Dirac delta function 
enforcing that we have always an eigenvalue pair (A, —A) of the Dirac operator if A 7^ 0. 
The two blocks are reminiscent of Vandermonde determinants and are artefacts of the 
zero modes at a = 0. 

The joint probability density of Z^w is 

p^{Z)d[Z]=c+{l + a^) 



a 

T 



V 



n 



1 + - K + 



\2 



X det 



l<i<n 



l<i<u 



n+u 



where 

gi~^\zi,Z2) 

gr{Xi,X2) 



gr{xi,X2)S{yi)S{y2) + gc{zi)S{xi - X2)5{yi + ^2) 



exp 



n 



2 n 2 , ^^&+ 



^^^{Xi+X2) + -{xi-xi) + 



9c{z) 



Xi - X2 

\xi - X2I 
li- — j- exp 



erf 



n{xx — X2) 



A7+ 



{xi + X2) 



n 



-X 



ny + 



nm6+ 



-X 



(2.14) 
(2.15) 

(2.16) 
(2.17) 
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^ (^_iy{u-l)/2+n{n-l)/2 (^\ (27r) ''/2n-'^'/2-n(n+.) "Q , "Q ^. , ^ (2.18) 

Notice that the one point weight (/(^^ of and of D^/ is apart from the sign of the 
hnear shift the same. Also the other distributions show similarities with each other. 

Comparing pw with we recognize the major difference is the determinant which 
replaces the Pfaffian. The reason is a broken permutation symmetry in the eigenvalues 
of -Dw- We have to symmetrize over the eigenvalues z^"^^ and z^^^ separately. Since 
the two point weight g2~^^ only couples z'^^'^ with z^^'' but not two eigenvalues of one 
and the same set the symmetrization yields a determinant. Another crucial difference 
of pw to p5 is the distinction of real and complex eigenvalues reflecting the non- 
Hermiticity of -Dw- Interestingly the complex conjugated pairs only enter the two 
point weight g2~^\ In the continuum limit the interaction of a pair of real eigenvalues, 
gr, is suppressed and the term for the complex eigenvalues, gc, enforces the pairing of 
non-zero eigenvalues, (^A, — «A), along the imaginary axis. Again a block resembling the 
Vandermonde determinant appears and is again a relict of the former zero modes. 

In the next two subsection we motivate the polynomials constructed in Sec. O For 
this we consider the fc-point correlation functions of and D^. 

2.1. The k-point correlation function of 

First, we consider the fermionic partition function of with A^f source terms, A = 
diag (Ai, . . . , AatJ, 

Zj;;"'-^(A) (X d[D^]P{D^)l[det{D, + \,l2n+u). (2.19) 

i=i 

The unit matrix of dimension 2n + z/ is denoted by l2n+i/- In the microscopic limit this 
partition function corresponds to the integral [201 EI] 

A \ n>l 

2^. 

U(Nf) 



C4 



^N^'^ ( tJt ) / d^{\J)d^\yU (2.20) 



X exp 



^ tr(f/ + f/-i) + i tr(A7]l^,. + A)(f/ - U'^) - ti{U^ + U'^) 



This is the effective Lagrangian of the Wilson-Dirac operator of the partition function 
with Af fermionic quarks with a degenerate quark mass mg and Af source terms AylArj+A 
proportional to 75, cf. Refs. [271 EHl EHl [30]. An integration over the variables friQ and 
A7 weighted by two additional Gaussian will yield the two low energy constants and 
Wj proportional to two squared trace terms (SDIEIIEII- Here we will not consider these 
integrals. 

Employing the joint probability density ps, see Eq. 02.61) . we combine the 
Vandermonde determinant and the characteristic polynomials to a quotient of two 
Vandermonde determinants. Then we rewrite the finite n partition function 02.191) 
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as 



(2.2i; 



2n+v 



X Pf 



\gi \x^,Xj)] {xl ^g[ 

^ l<i,j<2n+v l<i<2n+u 

l<j<u 



{-xr'g[-\x,)} 



l<i<i' 
\<j<2n+v 







We want to consider a little bit more than the partition function namely the /c-point 
correlation function. For this purpose we only integrate over 2n -\- v — k variables, 
X = diag {xk+i, . . . , X2n+v)- The remaining variables x' — diag {xi, . . . , Xk) are the k 
levels we look at, i.e. 



2n+v 



(x, -A) 



2n+v—k 



X Pf 



A.,(A) 

\g'f\xi,Xj)\ {xr^9^^\xi)}_. 



(2.22) 



{-x]r^gV{x,)} 



\<i<2n+v 
l<j<u 



l<i<v 
l<j<2n+v 







The idea is In the Vandermonde determinant of the numerator we can build an 
arbitrary basis of polynomials from order to order 2n -\- u -\- Ni — 1, 

A2.+.+^,(a;,-A) = (_l)(2n+.H-iV,)(2n+.+iV,-l)/2 



X det 



l<i<2n+v 



l<i<2n+iJ 
u<j<2n+u+N(-l 



l<i<Nf 
0<j<i^-l 



l<i<N( 
u<j<2n+i^+N{-l 



. (2.23) 



Also the entries of the Pfaffian can be transformed by adding rows and columns with 
each other, 

\gi~\xi,Xj)\ {xr^9[~\xi)}_. 



Pf 



l<i<2n+u 



l<i<v 
l<j<2n+v 







Pf 



{g[-\x,,x,)] {p^-\x,)g'i\x,)} 



\<i<2n-Vv 



0<i<t/-l 
\<j<2n+v 



(2.24) 



where we change the basis of the monomials to the polynomials ^ and the two-point 
weight 5f2~'' to G^^^ . 
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To shorten the notation we define the scalar product of two integrable functions /i 
and /2 with the one-point weight g['^^ 



(/l|/2)_(±) 
9l 



The same can be done for the two-point weight G2 
product 



(-) 



(/1I/2 



G. 



(-) 



det 



fi{z)f2{z)g^t\x)6{y)dxdy. (2.25) 

We define the anti-symmetric 

Gi'\zi, Z2)S{yi)5{y2)dxidyidx2dy2. (2.26) 



/l(^l) /2(^l) 
fl{z2) f2{z2) 



Both definitions are extended to the complex plane by Dirac delta functions because we 
want to discuss the situation for both random matrices and in a unifying way. 
In the next step we employ the de Bruijn-like integration theorem derived in 



Appendix A.l which yields 



<r^(x',A)cx 



-Pf 



-(m: 



^(12) 

(12) j^m) 



Ajv,(A) 

for the fc-point correlation function. The matrices in the Pfaffian determinant are 



(2.27) 



M 



(11) 



" {p'r'ip'f') 


0',-' 


W\ 








{it'ipy'] 


aV 


(«.'-•! 


'if'') 






. -(pT'lPf' 


\[-' 




wr 








(2.28) 







\x 




P\ \X3) 


pJ^Va,) " 


R 


\x 






r^(-A.) 




-p\ \xj)g[ ' 



















" 









































(2.29) 



(2.30) 



In the Pfaffian fl2.27p the indices i and j of the rows and columns are (0 . . . z/— 1, v . . . 2n+ 
z/ + A^f — l,0...z/ — 1... A^f) from top to bottom and left to right. Please 

notice that regardless what the polynomials and g^^) and the modified two-point 
weight 6*2" are Eq. fl2.27l) tells us that the joint probability density can also be written 
as a single Pfaffian. We have only to choose k = 2n + u to see that this statement is 
true. However the representation fl2.27p is quite cumbersome. A more compact one is 
given in subsection 14.11 

The aim is now to choose \ p[ ^ and G2 such that the matrix M^^^ becomes 
quasi- diagonal since we want to invert this matrix. A quasi- diagonal structure is 
equivalent to the conditions 

for 0<i,j<z/-l, (2.31) 
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H2i+u+l\H2j+u+l y (_) 

(-) I (-) 



0, for < 2 < z/ - 1, < j < 2n + A/'f - 1, 

0, for 0<z,j<i^-l, 

0, for < 2 < z/ - 1,0 < j < 2n + ATf - 1, 

0, for < z,j < 2n + iVf - 1, 

0, for < 2,j < 2n + A/'f - 1, 

o5"^5ij, for Q <i,j <2n + Ni-l. 
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(2.32) 
(2.33) 

(2.34) 

(2.35) 

(2.36) 

(2.37) 



The constants /ij ^ and o^- ^ are the normahzation constants of the polynomials. In 
Sec. [3] we will see that this system of equations have indeed a solution. We will give an 
explicit construction of them. 



2.2. The {kj- , k\) -point correlation function o/-Dw 

The next case we want to consider is the fermionic partition function of with N{ 
quark masses, m = diag (mi, . . . ,mN{), 



(X d[Dw]P{Dw)Y[det{Dw +mjl2n+u 



In the microscopic limit it corresponds to 



X exp 



\2nJ 



dfi{U)det''U 



(2.38) 



(2.39) 



U{Nf) 



i tr(m6l7Vf +m){U + U'') + y ti{U - U'^) - a" ii{U^ + U'^) 



This is the effective Lagrangian of the Wilson-Dirac operator of the partition function 
with iVf fermionic quarks with a non-degenerate quark masses meH-ATf + rh and one 
source term A7 proportional to 75 [271 EH [29l [30]. Again one can integrate over the two 
variables mg and A7 weighted by Gaussians to obtain the two low energy constants 



and but we will consider Z 



{n,v,+) 
Nf 



without these integrals. 



The partition function with the joint probability density pw reads 
d[Zl- 



oc 



. ^2n+u+N{{Z, —m) 



X det 



AiVf(m) 



(2.40) 



■' l<i<n 



{{xfr'gt\xfny^^)} 



l<i<v 



Since the permutation symmetry in the eigenvalues of is broken we have to 
consider a two parameter set of eigenvalue correlation functions. The number of 
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eigenvalues z^'^'> is independent of the number for the eigenvalues z^^\ Let Z' 
diag(zi'\ 



JO (1) 



^ ). Hence we define the {k,- , /ci )-point correlation function, 



p(n,j/,+) ^7/ 



/n n+v 
n n d^fdyf 



(2.41) 



j = kr +1 



i=fei+i 



X 7 — - — r clet 



•' l<i<n 



As in subsection 12.11 we construct an arbitrary basis of polynomials in the 
Vandermonde determinant to inflict some conditions on them later on, 

A2„+.+^,(^,-m) = (_l)(2n+.+A^,)(2n+.+7V,-l)/2 



X det 



p'*'(4") 



l<i<n 
0<j<u-l 

l<i<.n+i/ 
0<i<i'-l 

l<i<Afi 
0<i<i^-l 



l<j<n 
i^<i<2n+i/+Aff-l 

l<i<n+!^ 
l<i<Afi 



(2.42) 



Also the other determinant in the numerator can be transformed, 



det 



{#'(4", 4")} 



l<i<n 

l<j<n+i/ 



det 



{gW(..">,4")} „.< 

l<i<n 
J J 0<.i<.u—l 



The whole procedure works analogous to the one for D5, cf. subsection 12.11 
Let 

fM'^) /2(^(')) 



(/il/: 



2 M+) 



\ I det 



(2.43) 



be the anti-symmetric scalar product of two integrable functions /i and /2 with respect to 
the two-point weig ht Notice that G^^ as well as g^^ are not anti-symmetric under 

a permutation of their entries whereas the two-point weight G2 ^ is anti-symmetric. The 
reason for this is again the breaking of the permutation symmetry in the eigenvalues of 

Considering the {k^ , k\)-point correlation function we apply the de Bruijn-like 
integration theorem derived in Appendix A. 2 to the partition function fl2.40p and find 



R 



Pf 



M 



(11) 

(12) 



M 



(12) 

+ 

r(22) 



(2.45) 
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) 




























(2.46) 



(M 



(12)nT 






































(22) 













" 


























































. 



(2.47) 



(2.48) 



with = dxdy. In the Pfaffian 02.451) the indices i and j are in the range 

(0...z/-l,z/...2n + i/ + A^f-l,0...i/-l,l...fcr,l...fci,l...fcr,l...fci,l...A/'f) from 
top to bottom and from left to right. 

Please notice the similarity of Eq. fl2.45p with Eq. fl2.27p . \i = n and ki = n + u 
the correlation function is equal to the joint probability density pw . Hence pw can also 
be written as a single Pfaffian which can be cast into a compact form, see subsection 14.21 

As in subsection 12.11 we want to invert and, thus, quasi-diagonalize the matrix 
M^^\ This yields the following system of equations 



G. 



{+) 



(+)| J+) 



Pi h 



G\ 



(+) 



{+) I (+) 



G. 



{+) 



(+) l^(+) 



4+1.1^2 



hf^bij, for 0<i,j<z/-l, 

0, for < 2 < z/- 1,0 < J < 2n + iVf - 1, 

0, for 0<2,j<z/-l, 

0, for < z < z/ - 1, < J < 2n + iVf - 1, 

0, for < i,j < 2n + A^f - 1, 

0, for 0<i,j <2n + A'f-l, 

ot^bi^, for < i, j < 2n + A'f - 1 



(2.49) 
(2.50) 
(2.51) 

(2.52) 

(2.53) 

(2.54) 

(2.55) 



with the normalization constants /i^'*''' and o^' ■ Comparing this system of equations 
with the one of we recognize that they are of the same form. Hence, if we solve them 
in a general setting we solve them for both random matrices and • 



(+) 
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3. Construction of the polynomials and some of their properties 

In subsection 13 . II we construct the orthogonal polynomials as well as the skew-orthogonal 
ones starting from the conditions fl2 .31112.371) and (12.49112. 53]) . Furthermore we give 
explicit expressions of the modified two point weights, G'^^ and specify the constants hi 
and ol^2i- Recursion relations of the polynomials are shown in subsection l3.2[ With help 
of these relations we derive the Christoffel Darboux-like formula, in subsection 13.31 In 
subsection [331 we rewrite the polynomials as well as the Christoffel Darboux-like formula 
to random matrix averages and take the microscopic limit of them in subsection 13.51 



3.1. The polynomials 

The starting point of the construction are the monomials 



xnj[z) 



with j G Nq. 



(3.1) 



With help of the general formula for the orthogonal polynomials of the one-point weight 
(^l^^ |44j as a quotient of determinants we find 

|(mi|mj)^(±) 



det" 



(mi|mj)^(±) 



det 



0<i<l-l 
0<j<l 



0<j<l 



(3.2) 



in monic normalization, i.e. p] 



z'' + . . . Since gi is a shifted Gaussian in our case 



the polynomials are shifted Hermite polynomials Hi in monic normalization. 



Pi 

P?\z) 



a 



2\ in 



n 



- Hl\J-zT\-^'l 



(3.3) 



This agrees with Refs. [20| [2T| [22| [23] where a mixing of the eigenvalue statistics with 
a finite dimensional CUE was found. The normalization constant is 



a 



2 \ '+1/2 



^2n { — ] n exp 



2n 



(3.4) 



Thus the normalization constants of the orthogonal polynomials iare the same for D-w 
and D^. 



Starting from the orthogonal polynomials we want to construct the polynomials g. 

,(+) 



(-) 

u+l 



and ql_^_\ fulfilling the orthogonality conditions (12.321) and fl2.50p . respectively. Let A^f 
be even for simplicity. If it is odd the antisymmetric matrices M^^^ are never invertible 
because their dimensions are odd. In such a case we extend the partition function by 
one fermionic fiavour and remove it at the end of the day. 

As for the orthogonal polynomials we begin with an intuitive definition 



X Pf 



Pf 



(3.5) 



i'<i,j<u+2l-l 



92 



Pi^\^) 



v<i<v+2l 



-pf\z) 
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-1 



I J±) 
92 



(3.6) 



xPf 



^ ^2^^' j i'<i,j<i^+2«-l 



/ „'-^> I ;/<i<iy+2t — 1 ;/<j<i^+2t — 1 



52 



'Pi+2l+li^) 



Pu+2l+li^) 



-pfi^)\ -p)^,'.,,Az) 

J i/<j<iy+2«-l 

which is similar to the ansatz of the skew-orthogonal polynomials in xRMT with non- 
zero chemical potential and Dyson index /3 = 1,4, see Ref. The ant i- symmetric 
products of g2 and g^2'^ are defined similar to Eqs. f l2.26p and f l2.44p . respectively. 
One can readily prove that the orthogonality conditions fl2.32p and fl2.50p are fulfilled. 
The multilinearity of the Pfaffian allows us to pull the scalar product into the Pfaffian. 
Then one row and one column is zero. For example, let / G {0, 1, . . . , z/ — 1}. Then the 
orthogonality of the polynomials p^^'' yields 



Pf 



\Pl \'iu+2b/g[^^ 



9i 



^<i<v+2b 







Pf 



9-2 



Pf 



p. 



{0} 



v<i,j<u+2b-l 

{0} 

u<i<v+2b 




Pf 



Pf'^\pf\ (±) 



92 J v<i,j<v+2b-l 



0. 



(3.7) 



In a similar way one can prove the other relations. The normalization constants in 
Eqs. f l3.5p and f l3.6p are finite since they are proportional to the quenched partition 
functions (ICTj) and f l238|) . i.e. 



Pf 



oc Z, 



^0 



(3.^ 



v<i,i<v+2l-l 



It can be easily shown that the polynomials gf ■* are in monic normalization, too. The 
constants f l3.8p are related to the ones in Eqs. 02.371) and 02.551) by 



Pf 



' >9'^\ 



y<i,j<v+2l+l 



Pf 



(3.9) 



v<i,j<y+2l-l 

Combining this identity with relation f l3.8p the constants of'' are mostly the quotient 
of two quenched partition functions. 
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The constants O;^"* can be directly calculated by the normalization of the joint 
probability densities (12. 6p and f l2.14p . In particular we have the two identities 

u—l n—1 

{-in2n + u)\l[h,l[o^-^ 

j=Q j=0 



(3.10) 



1 



(1 



a 

Y 



2 , + ^ 2 



n 



mi 



and 



1 



^).(.-l)/2+n(n+l)/2^,(^ + Z/)! J] J] of^ 

j=0 j=0 
2 



(3.11) 



;i + a2)"(«+-i/2)an+-%xp 



a 



2 n + z/ 2 
/i, + /ii 



X2 



These identities can be derived by using the de Bruijn-like intergation theorems in 



Appendix A With help of Eqs. (13.101) and (13. lip we conclude 

2\ 21+u+l 



-411(1 + iy)\ 



TT 



n(l±a2) 



n 



a exp 



2n 



it^r +^i) ^(/^r Tt^i] 



(3.12) 



Hence the normalization constant is linear in a for small lattice spacing and is 
proportional to n~'^^~'^~'^l\{l + z/)! in the microscopic limit. 

The polynomials q~^2b ^"^^ ^^u+2b+i ^^^^ orthogonal to pf''' , / G {z/, . . . , z/ + 26} 
and / G {z/, ...,z/ + 26— l,z/ + 26+l}, respectively, corresponding to the two-point weight 
g^'^ since the rows and columns are not linearly independent anymore. For example 



Pf 



„(±)|j±) ^ 



I u<i<u+2b 





Pf 



= 0, (3.13) 
the Ith row and column and the last ones are the same. In the same way one can prove 



the skew-orthogonality of p\^^ , / G {u, . . . , u + 2b — 1, u + 2b + 1} , with 
to the definitions (13. 5p and (13. 6 p the polynomials g^^^ are a linear combination of p^^^i 
with < / < 2n + A^f — 1. Therefore the polynomials are indeed skew-orthogonal with 
respect to g^2^\ In particular they fulfill the conditions similar to the relations (12.351 



(±) 

v+26+1- 



Due 



I2.37P and (12. 53112. 55|) by exchanging 6*2^ "* — )■ g^' ■ However the remaining conditions are 
not fulfilled. This is the reason for modifying the two-point weights. 

The simplest way to enforce the remaining conditions is the projection of the 



measures (y'g^'' onto the polynomials q]^'i only. This means the polynomials p\ 
< / < z/ — 1, have to be in the kernel of G^^\ We make the ansatz 



(±) 



{flj2)(j(±) 



ifl, f2)j±) 
92 
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+ E (3.14) 

i,j=0 ^ * 

Indeed we have 

[p\^\ f) (_) = (p!^^ /)g(+' =0' functions / and < / < z/ - 1, (3.15) 

(#^.^%,^(#^.n,.. (3.ia) 

Thus all orthogonality conditions are fulfilled. 

The explicit expressions of Gi~^ and Gi~^^ are 

Gi'\xi,X2) = gi~\xi,X2) (3.17) 



- E/^ /p5"^(^')^^^(^>2)rfx'p,(xo(-)^?!-^(xi) 

+ [ p\~\x')gi-\x,,x')dx'p^-\x2)g[~\x2) \ + J2 

X j p!"^(a;i)p5~^(4)^2~^(a;i,4)^^[a;V!~^(^i)Pi"^(a;2)^l"Ha;i)^!^^(x2), 
G^+)(^i,Z2) =gt\zi.z2) (3.18) 

- E/T (^75^^(^i,.')-^?^^H^',^i))^[^>r^(^2)^7!^^(^2) 



+ jp^\^>f\^'2){9t\^[.^'2)-9t\^'2.A))d[z'] 

Xpr^(;^l)p5+^(^2)^?S+^(^l)^7^H^2). 

The change of the two-point measures is restricted by linear combinations with other 
rows and columns in the Pfaffian fl2.24p and the determinant (12.421) . Essentially we 
add the orthogonal polynomials pl^'' to the weight. Thereby we have to recall that 
everything which is done with the rows has to be done with the columns in the Pfaffian. 
This is the reason why G2~'' stays anti-symmetric whereas (^2'*'^ is asymmetric in the 
entries. 

3.2. Recursion relations 

The recursion relations of the orthogonal polynomials are 

-^{x) = Ip'tlix), (3.19) 

2 1 2 

xpl^^ (x) = (x) ± (x) + — (^) • (3-20) 

1% Th 
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They result from the orthogonality relation fl2.3ip and the two identities 

(/^^^Vi|/2),(±) = -(/i|/^(^)/2),(±), (3.21) 

(mi/i|/2)r±) = (/ilmi/2) r±, ' (3.22) 

for two arbitrary integrable functions /i and f2- The function mi is the monomial of 

order one and the differential operator D'^^'> is the creation operator of the harmonic 
oscillator corresponding to the measure g^\ 

Identity f l3.22p cannot be extended to the measures G^'^ and G'^'^ or equivalently g2~^ 
and ^'2'''^ i-e. 

(mi/il/2) r±, 7^(/l|mi/2) r±). (3.24) 

92 J2 

However Eq. (13.211) has an analogue. Defining the differential operator 
one can readily verify 

(5'*'/.IA)^,„ = - (/.|5<*'A)^.„ . (3.26) 

The starting point of such a proof is the differential equation 



d n 



gt\zi,Z2) = Q. (3.27) 



Notice that the differential operator is restricted to the real eigenvalues of and 
and to the real part of the complex conjugated pair of 2; = = ^2 of due to the 
Dirac delta-functions. 

Both operators D'^^^ and D^^^ are closely related with each other which is quite 
advantageous. For example the action of Z)*^^^ in the scalar product fl2.25p is 



{D^^^hm^m = - + ^^4^^(/i|/2),(±), (3.28) 



/ir =F /il 

Jl\J2)g[±) + 

where fl^ is the first derivative of /2. 

We consider the action of D^^^ on the polynomials p^{^^ and q'^i- The recursion 
relations flCT]) and f lX^ yield 

5<±)pr>(.) = - + f^^p-fiz). (3.29) 

The polynomial D^^^g^^) can be expanded in the polynomials {pj" 

5(^)£i(.) = i: «r gS(.) + E (.), (3.30) 

where a)^'' and are the coefficients which have to be found. 
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In Appendix B we derive the recursions 



(±) ( 



(3.31) 

(±) 

^^itL2(^), (3.32) 
o 



a 



/-I 



with the coefficients 



~(±) 



(2/ + l)^^. 



2 



(/ + l)2(/i, T/ii)2a2 



a 



^!^ + 2« + l hy+2l-l 



(3.33) 
(3.34) 



The recursion formula for q)^j^2i+ii see Eq. (13. 32 p . is restricted to / > 1. For / = we 
have to omit the last term, i.e. we set g^^2 = 0- This formula is quite useful to find 
Christoffel Darboux-like formulas, see subsection 13.31 



3.3. A Christoffel Darboux-like formula 

For the calculation of spectral correlations the Christoffel Darboux formula is quite 
useful. However searching for such a formula of skew-orthogonal polynomials proved as 
a difficult task [64J. The same is true for the polynomials g^^] for which we want to 
simplify the sum 



n-l 



SS(^1,^2) = Xl^^^* 



/=0 



qy+2ii^'i) qu+21+ii^'i) 



(±) 



(3.35) 



For the orthogonal polynomials p^{^^ we already know such a result 



^ 1 



- (i)/ X (±)/ X 1 P^^^^ {Zi )pS (z2 ) - Vv"^ {Z2 )pS {Zi ] 

«=0 



(3.36) 



Zl - Z2 

Identity fl3.36p is a direct consequence of the three term recursion relation fl3.20p . Hence 



we pursue the same way for Eq. fl3.35p which is done in Appendix C The Christoffel 
Darboux-like formula for the skew-orthogonal polynomials is 



'^n-\{Zi, Z2) 



n 



a^o. 



(±) 

n-l 



exp 



n 



n 



[Zl + Z2, 



X det 



£2n-2(^l + £L(^1 + 

?S2n-2 (2:2 + X) gj,+^2n (^2 + x) 



dx. 



X 



(3.37) 



This result only depends on a few polynomials as it is already well known for the original 
Christoffel Darboux formula, cf. Eq. f l3.36p . 



3.4. Representation as random matrix averages 

As we have already seen in subsection 13.21 all skew-orthogonal polynomials are easy 
to derive if we know a compact expression for I even. For this purpose we want to derive 
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a representation as an integral over a random matrix. For the orthogonal polynomials 
the well known expression of this kind is 



n 



2\ I 



j (lei{zti - H)exp 



n 



2a2 



tr if ^ — /ii 



n 



d[H], 



(3.38) 



where H is a I x I Hermitian random matrix with the measure 

d[H] = Y[dHii Y[ 2dRe Hijdlm Hij. (3.39) 

i=l l<i<j<l 

This random matrix integral can be drastically reduced to a small number of integration 
variables by the supersymmetry method [511 ESI [SH [56] . A famous representation of the 
Hermite polynomials can be calculated in this way, 

27r 



Pi = ^ / ^^p 



In \ n 



(3.40) 



The corresponding Rodrigues-formula for the Hermite polynomials is a simple lemma 
from this, i.e. 



a 

n ^ 
2\ V2 



V 


■ J2 




1 exp 







n 



2a2 



(3.41) 



Hi 



The aim is to find the formulas analogous to Eqs. (13. 38113. 4T1) for q'j^2i- 

We compare the definition (13. 5p with Eqs. (I2.27P and (12.451) for = fcr = k\ =0, 
n = I and iVf = 1. Since we were free of chosing the two point weight G^2^^ and the 
polynomials at this step of the calculation, Eqs. (I2.27P and (I2.45P are also valid 
when replacing by g^2^^ and by p^^i- Moreover, q^^2i ^^^^ equal to 



(_l)'^('^+l)/2 



Pf 



v-1 

n h, 



xPf 



Pflpf) 

0<i<i^+2« 

{-P,(^)W} 



(3.42) 



9^*' J 0<i,j<u+2l 



(±) , 

N-f-* I^J /gi o<i<i.+2« 



0<j<!/+2i 



o<i<i'-i 



0<i<i'+2i 













Then the polynomials q^2i partition functions with one fermionic flavour 



(3.43) 
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By means of the supersymmetry technique [5T1 [53l [SU [56] we derive the result 



(±l)'+^Z!(/ + z/)! 



X exp 



n 







exp 


( 


[0,27r]2 










V n 



Tz e 



n 



(3.44) 



m 



Appendix D Notice the similarity with Eq. fl3.40p . 



Again we can ask for a Rodrigues formula and indeed it is a direct consequence of 
Eq. ^M)- We find 



[±l)i+^ I — 
n 



a 



2 \ 



Ql 

83 



exp 



x=0 



a 

2n 



exp 



X + 



n 



(3.45) 



n 



n 




X 



Performing the derivatives we find an explicit expression in terms of Hermite polynomials 
for the skew orthogonal polynomials g^^^' 



= (±1) 



l+v 



1 



X H, 



u+j 



± 



n 



n 



3=0 



l\{l + v)\ 



j!(Z-j)!(^ + j)! 



i2i 



(3.46) 



z — 



z — 



The polynomials for ql^2i+i t)e readily obtained with help of relation f l3.3ip . 
Remarkably the prefactor of the single summands is exactly the same as the one of 
the modified Laguerre polynomial lJ^'' when replacing the Hermite polynomials by 
monomials. 

The limit a — )■ yields the generalized Laguerre polynomials, 



a->-0 



1 

n 



n ' 



(3.47) 



n 



±- 



n 



which is in agreement with Ref. The large a limit with fixed variables ^Jn/a?z and 
^faFJuy^x /r is a product of two Hermite polynomials 

9N (2« + l/)/2 



a>l / 



n 



n 



v+l 




n 




(3.48) 

Both limits can already be directly derived from the random matrix model, cf. Eqs. (12. ip 
and (12. 2p . For a = we have a xGUE whose orthogonal polynomials are the Laguerre 
polynomials. Recently it was shown that the xGUE has also a non-trivial Pfaffian 
factorization whose skew-orthogonal polynomials of even order are the orthogonal 
polynomials itself. Hence the limit (I3.47P agrees with the observation in Ref. [62] . 
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In the large a limit the off-diagonal blocks W and IV^, see (12.11) . are suppressed. 
Therefore we end up with two decoupled GUE's. One is of dimension I and the other 
one of dimension I + u. This indeed yields a product of two Hermite polynomials, cf. 
Eq. ([S3S1). 

A particular case of the polynomials can be obtained for the random matrix D^. 
Let = —f'l = fJ' and a = 1. Then we have a. 2n + u dimensional GUE, cf. Eq. (12. 2p . 
Indeed we also get the corresponding Hermite polynomials. Equation (I3.44p simphfies 
to 



{±iy+n\{i + u)\ 



X exp 



n 





[0,27r]2 


--)( 




n/ 






I \/nz 



exp 



2n 



(3.49) 



for the even polynomials and 



\- Z 

n 



(3.50) 



\/nz -= 

n 



for the odd ones. Therefore all polynomials are given by Hermite polynomials 
corresponding to the same Gaussian distribution. 

Another useful random matrix integral representation would be the one for the 
Christoffel Darboux-like formula (I3.35p . For the orthogonal polynomials pf""* such a 
representation is well known, 



Z\ - Z2 



n 



2 \ ('^-i) 



det(2;il,^_i — H) det(z2lt/-i — H) exp 



2 \ 2' 

tr ( if ^ —^1 



2a 



n 



(3.51) 
d[H] 



with a (z/ — 1) X (z/ — 1) Hermitian matrix H . With the supersymmetry method 
[STj [531 [SU [56] one can also find the representation 

V'^\z.)pf\{^2)~pf\z2)vf\{^.) 



Z\ - Z2 



(3.52) 



v\{v-\\\ 



exp 



2n 



-— tr f/' + tr diag {z^, z^) T -^^2 U 



n 



det 



-v+l 



Udfi{U), 



U{2) 



where dfi{U) is the normalized Haar measure of the unitary group U (2). 

In Appendix E we show that the Christoffel Darboux-like formula (I3.35P is 
essentially the partition function with two fermionic flavours, i.e. 



Zl,Z2 



[Zl - Z2 



-Zl, -Z2 



(±) r^{n,u,±) 
On Z/n 



(3.53) 



^n ^0 
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Also the two flavour partition function can be mapped to an integral over unitary 
groups by performing the same calculation as for the one flavour partition function, see 
the discussion above Eq. (13.441) . Therefore T^lt^ is an integral over a compact set, 

21+u+l 



eW(zi,Z2)= - 



(n + l)\(n + u + 1)\ /n(l±a2 



TT 



n 



1 ±a2 



1 

a 



X exp 



/ 2 2\ it / 



2n 



X / diJ,{Ur)diJ,{Ui) exp 

U(2)xU(2) 



- — (tr + tr U?) + -trUr Ui 
2n n 



X exp 



X exp 



tr ( —^^^2 - diag (zi, Z2) ] 
tr ( ^-|^l2 T diag (zi, 2:2) ) Ui 



det-"f/. 



(3.54) 



Equations (I3.44p and (I3.54p are suitable for discussing the asymptotic behavior as it is 
done in subsection 13. 5[ 



3.5. Asymptotics 



The microscopic limit (n — t- 00, see dissension after Eq. (12.111) ) directly relates chiral 
random matrix theory with QCD in the e-regime. Hence we want to know the 
expressions of the polynomials as well as the one of the Christoffel Darboux-like 
formula S^'j'i in this limit. 

For an arbitrary function / which is n-independent and smooth on the group 
U (fc) X U (k) the following asymptotic result exists 



f{Ur , Ui ) exp [n tr U,Ui - n tr MJ^ Ui] d^i{U, )d^{Ui ) 



(3.55) 



U(fc)xU{fc) 



fe-1 „ 

(2vr)-^/2n-^'/2e"'^n^- j fiU,U-')df,{U). 

■^'^^ U (fc) 



Equations (I3.44p and (I3.54p are particular cases of this identity. Hence we have 



(±) / _Z_ \ n»l 



-l)^(±l)"+^0^e-" 



exp [-a^e'^-^ + e 



'i2ip 



(3.56) 



[0,27r] 



X exp 
for the polynomials and 



1 



^(±) I ^ f2_ \ n^l 

2n2nJ %^/2ti a 



exp 



- ^ T 4a^ 



[Zl - Z2\ 



(3.57) 



X j exp 

U(2) 



a'(tr U' + tr U'') + - tr (mg + \7)t2 - diag (^i, $2) U/ 
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X exp 



tr (^{fhe - X7)l2 =F diag (zi, U ^ 



defUdfi{U). 



for the Christoffel Darboux-like formula, cf. Eqs. f l2.20p and (12.391) . In both equation 
we apphed Stirhng's formula to the factorials. 

In the case of the polynomials q^2n{^) '^^ ^^^^ able to integrate over the domain, 

^''^'^^^}!Zf''^~^ f; / d^id^^ew [-a^(e- + e"-)] (3.58) 



(27r)3/2 



j=-oo. 



[0,27r]2 



X exp 



me + A7 - z) e*'^2 + - (me - Xj Tz) e 



gl!y(/P2g«j(vi-2V2) 



rhe- XjTz 



u/2+j 



j=-oo 



TUq + X7 - Z ^ 

X I„+2j { \l {ma + A7 - z){me - A7 =F S') ) Ij{-a^). 



Due to the modified Bessel functions of the second kind Ii{z) oc (2;e/2|/|)l'l/ A/27r|/| oc 
(z/2)l'l/|/|!, for |/| ^ 1, the series rapidly converges and is numerically more stable than 
the integral (13.561) in simulations. 

Unfortunately it is much harder to find such a sum for the Christoffel Darboux-like 
formula. However we can diagonalize the unitary matrix U and find 



v(+) / ^ ^ 1 ^ 



exp 



A? 



d(fidip2 sin^ 



(fl - (p2 



(3.59) 



[0,2tt]' 



X exp 



.i=i 



2a cos ipj + 



4a 



+ ivipj — iXj sin Lfj 



X 



exp \zi cos v?! + % cos ^92] — exp \z2 cos + zi cos 932] 

cos — COS(y92 



for and 

v(") ( ^ ^\ 
'2n'2n 



4(27r)5/2a 



exp 



dipid(f2 sin^ 



(3.60) 



[0,27r]: 



I I 2a sin — | + — ifriQ cos 



X exp 

exp [I'zi sin + iz2 sin (y92j — exp [i'z2 sin y^i + iz'i sin ^92] 

X 

sin Lfi — sin (y92 

for D5. These two formulas are quite suitable for the applications discussed in Sec. IH 



4. Application to Wilson RMT 



The results of the previous sections are helpful to simplify the /c-point functions of as 
well as of D^J^/ . A Pfaffian factorization of the eigenvalue correlations of was already 
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given in Ref. [23]. We obtain this structure in Sec. 14. ![ too. Moreover we express the 
kernels of the Pfaffian in terms of two flavour partition functions which has proven 
fruitful in other random matrix ensembles, see Ref. [57] and the references therein. The 
unquenched (fcr, ^i) -point correlation function of is shown in Sec. 14.21 which is a 
completely new result. Also this result displays a Pfaffian factorization whose entries 
are two flavour partition functions. 



4.1. The Hermitian Wilson random matrix ensemble 

In the /c-point correlation function fl2.27p we encounter an integral transform of the 
orthogonal and skew-orthogonal polynomials, cf. Eq. f l2.29|) . Thus we define the integral 
transform of the skew-orthogonal polynomials. 



(iu+i{^)^2 \x,x)dx 



(4.1) 



(ll+i{^)92 '{x,x)dx-2_^ p) '{x)g\ '{x) 



j=0 



hj 



li+ii^)^2 \x,x)dx. 



The same integral transform for the orthogonal polynomials p\ \ < I < u—l, vanishes. 



(4.2) 



cf. Eq. fl3T7D . 

Using the identity 

A B 
-B^ C 



Pf 



Pf APf [C + B'^A-^B], 



(4.3) 



where B and C are arbitrary and A is invertible, the A;-point correlation function with 
an even number of fermionic flavours A^f = 2nf, see Eq. fl2.27p . is 

\A;(fc+l)/2 

(4.4) 



Pf [A'g ' ^ '^''( — Aj, — Aj)]l<4j<2nf 







-,n+nf) 


Xj, Xj) 


-At-"' 


\Xj , X j ) 


-;r(-'"+"')(-A„a;.) 


X Pf 




-,n+nf) 


Xj, Xj) 


At' 


^Xj, Xj) 


i^r-+-^(x.,-A,) 




At 


n+rif) ^ 


Aj , Xj ) 


aS->( 


Aj, Xj) 


-A,) 



which is the main result for the Hermitian Wilson random matrix D^. The indices i 
and j of the Pfaffian in the denominator take the values (1, . . . , fc, 1, . . . , A;, 1, . . . , 2ni). 
The functions in the entries are 



ir^'"+"'')(Xi,X2) =G5"^(Xi,X2)+ iq^^* 



i=0 



tUi(^2) gi;L(^2) 



(4.5) 
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u-l 



(— ,n+nf) 



{Xi,X2] 



= E 

1=0 
n+rif — 1 

+ E 



1 

hi 



K. 



(-,ra+nf) 



Xi,X2] 



n+Tif — l 

E 

z=o 



1 



1 



det 



det 



{-) 



Xi) 



(-) 

'iu+21+1 
%+21+li^'i 



(lu'+2li.^2) 



(4.6) 



= - Si+lf-i(a;i,a;2). (4.7) 

The fc-point correlation function for an uneven number of flavours can be derived by 
shifting one of the variables A to inflnity. Then we get in one row and one column 
of the numerator and the denominator of Eq. 04.41) the skew-orthogonal polynomial 
ql+2{n+n,-i) and its integral transform qiX^^^^^_^y 

The case = 1 is the normalization. For k = 2n + u we have a compact 
representation of the joint probability density as a single Pfaffian determinant. 

The representation (14.5114. 71) in terms of the Hermite polynomials and the skew- 
orthogonal polynomials g^^) can be easily interpreted. The u former zero modes are 
broadened by a GUE of dimension u. The skew orthogonal polynomials can be identifled 
with the remaining modes and describe the spectral density thereof. Both spectra, the 
one of the GUE and the one of the remaining modes, are coupled by the sum in Eq. (14. ip . 
They manifest the repulsion of the former zero modes with the remaining modes which 
is given by the Vandermonde determinant in the joint probability density (12. 6p . 



The kernels 



(-,71+rif) 



and K. 



(-,n+nf) 



can be expressed in partition functions, too. 



In Appendix F we derive the following results 
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TT X2 - Xis-^O \ det(L'5 - {X2 + ie)l2n+u+2) / 
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f{x)dx f ef{x)dx 



Im 

£-i>o / X — le 



X^ + 6^ 



and 



(F(Dw))^..= (i^(75/^5)) 



N,u 



F{D^)P{D^)d[D 



w 



(4.9) 



(4.10) 



(4.11) 



for two arbitrary sufficiently integrable functions / and F and the definition of the 
probability density P in Eq. (12. 2p . The random matrix on the right hand side of 
Eq. (14. lip has the dimension {2N + u) x {2N + u) with index u. Hence we have to 
take the averages (14. 8 p and (14. 9 p over a Wilson random matrix with N = n + Uf. 
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Considering Eqs. fl3.53p . (14. 7p . (14 .Sp and f l4.9p we traced the unquenched A;-point 
correlation of back to partition functions with two fermionic, two bosonic and 
one fermionic and one bosonic determinant. Hence the structure of the eigenvalue 
correlations is in the same class of matrix ensembles as the /3 = 1 and /3 = 4 
standard ensembles, e.g. GOE, GSE, the real and quarternion Ginibre ensemble, the 
chiral GOE and the chiral GSE, see Ref. [57] and the references therein. When taking 
the continuum limit, a — )■ 0, the Pfaffian determinant will persist though we have then 
chiral GUE. This observation agrees with the result found in Ref. [62]. Therein a non- 
trivial Pfaffian was derived for all random matrix ensembles corresponding to orthogonal 
polynomials. Exactly this structure carries over to the finite lattice spacing result (14. 4p . 

The kernel K2~'"'\x,x) is equal to the quenched one point function of D^, denoted 
by P5(x) in Refs. [201 EI]- Due to the prefactor l/(xi — X2), see Eq. (14. 9p . we have to 
apply I'Hospital's rule which exactly agrees with the common definition of p^. 

The Pfaffian factorization (14.41) was already discovered in Ref. [23] but we made 
the connection to two fiavour partition functions. Furthermore the structure as well 
as the expression in two fiavour partition functions carry over to the microscopic limit. 
In this limit Wilson random matrix theory is directly related to the e-regime of Wilson 
fermions in lattice QCD [271 [2H1 [291 EOl [39] . Hence we found a neat representation which 
drastically simplifies the numerical realization of the /c-point correlation functions. 

The microscopic limit of the kernel '"^ is shown in Sec. 13.51 see Eq. (I3.60p . A 
derivation of the other kernels as well as a qualitative discussion of the results will be 
done elsewhere [63] . 

In the notation of Refs. [201 1211 122] the kernels are proportional to the two fiavour 
partition functions of the chiral Lagrangian 

k[~'°°\xi,X2) oc (xi - X2) Im Zo/2(m6,m6; A7 - xi, A7 - X2;as,aQ/j = 0), (4.12) 



in microscopic limit. Please recall that x = 2nx is fixed. The constants are essentially 
the product of the lattice spacing a times the square roots of the low energy constants. 



\/Wi, [201 1211 122]. We get the case ag/y 7^ when we multiply the expression (14. 4p with 

the partition function of A^f fermionic fiavours cancelling the Pfaffian in the denominator. 
Then we have to integrate over Gaussian distributions of mg and A7. Finally we divide 
the result by the partition function of Nf fermionic fiavours with Oq/j ^ which is also 
the two Gaussian integrals over mg and A7 of the partition function with a^/j = 0, cf. 
Ref. [M]. Please notice that we will lose the Pfaffian factorization when going from 

^6/7 = to 06/7 7^ 0. 



£l-i>0 

e2-s-0 
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(4.13) 



'°^'{xi,X2) oc (Xi - X2)^2/o("^6,"^6; A7 - Xi, A7 - X2; 08,06/7 = 0) 



(4.14) 
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As in the Hermitian version we define the integral transform of the skew-orthogonal 
polynomials Q'f,^). However we have to distinguish between left and right transformation 
because G^2^'' is not anti-symmetric anymore, 



u+l 



{I}gi'-\l,z)dm-J2 

j=0 

iI)Gi-'\z,I)dm 
{I)gi-'\z,I)d[^. 



(4.15) 



(4.16) 



Another difference to the Hermitian case is a non- vanishing integral transform of the 
orthogonal polynomials 

p[^\l)Gi-'\z,z)d[I\ = I p[^\l)Gi-'\z,m^ 
c c 
for < / < z/ - 1 due to Eq. (KWi . 

Again we consider an even number of fermionic flavours. Then we arrive 
main result for the non-Hermitian Wilson random matrix which is the [k^ , k\ 
correlation function f l2.45p . 
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(4.19) 
(4.20) 
(4.21) 
(4.22) 
(4.23) 
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where the indices i and j take the values (1, 
right and top to bottom. The functions are given by 



,k\,l,. . . , 2nf) from left to 
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(4.27) 



(4.28) 



(4.29) 



The result for an odd number of flavours can again be obtained by taking the limit 
of one mass to infinity. This will yield that one row and one column only depend on 
the skew-orthogonal polynomial ?^+2(n+nf-i) integral transforms 5l+2(n+nf-i) 

For kj- = fci = we find the normalization and in the case {k^,ki) = (n, n + u) 
wc have a representation of joint probability density ^ Pfaffian similar to the 

one of P5. Additionally, we can consider the particular cases {kr,ki) = (n, 0) and 
{kr ,k\) = (0, n + 1/) which are the joint probability densities for the eigenvalues z^^^ and 
z^'^ separately. These two joint probability densities are the ones for the right handed 
and the half of the complex modes {z^'^^) and for the left handed and the other half of 
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the complex modes (z^^^)- 

Again we recognize a natural splitting of the spectral properties. There are those 
terms, the sums with Hermite polynomials p['^\ which describe the broadening of the 
former zero modes. They are again manifested by the same GUE which we found when 
discussing and are located on the real axis only, notice the Dirac delta functions. 
Moreover we have the terms for the remaining modes given by the skew-orthogonal 
polynomials, g^^). The corresponding eigenvalues to these modes do not necessarily lie 
on the real axis. On the contrary the most eigenvalues are distributed in the complex 
plane, see Ref. [igiMlIM]. 

There is an interaction between these two kinds of spectra in the integral transform, 
cf. Eqs. f l4.15p and f l4.17p . This interaction directly follows from the Vandermonde 
determinant in the joint probability density fl2.14p . The repulsion obtained by this 
coupling effects the spectrum located on the real axis as well as the complex one. 

As for Dr, all kernels of the result f l4.18p can be traced back to two flavour partition 
functions. For the kernels Eq. (^2M, 7^:^^+'"+"^^ Eq. flCTjl . and 

Eq. fl4.29p . we know already appropriate expressions. In Appendix F we derive the 
results for the other kernels, 

Air(+'"+"')(zi, Z2) = Z2) - z,) (4.30) 

— 92 1^2, Zi) + 
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e^O \ det(Z;w - Xil2n+^+2 " ^^Ts) / n+nf,u 

Airl+'"+"*)(zi, Z2) = irl+'"+"*-)(^i, z,) - Air(+'"+"')(zi, z;) + AKf (^2, zl) 

+ ir;+'"+"'^(z*,4) (4.31) 



X 



(xi - X2)6{yi)6{y2) Im 

ei-!>0 
1 



n+nf,u 



^det{Dy^ - xil2n+v+2 - ^£llb) det(L)w - a;2l2n+!.+2 - ^£^275) , 
Air(+-"+"')(zi, z,) = Kt'''^-^\z,, z,) - K^^^-^''\z,, zl) (4.32) 
= l^Im / det(Z^w-.iW2) \ 

T\ Z1-X2 \ det(L'w - X2^2n+v+2 " ^£^75) / n+n,,v 

The kernel Ai^4^'"^"'^'' describes the correlation of the chiral distribution over the real 
eigenvalues with itself. This can be seen by the 75 weight of the imaginary increments 
and the Dirac delta functions, cf. Refs. [22]. The other two kernels fl4.30p and (14.321) 
represent the interaction of the chiral distribution over the real eigenvalues with the 
remaining spectrum describing the additional real modes and the complex ones. 

The quenched one point functions presented in Refs. [221 [13 Ell are given by 
the kernels k'^^'"^^ {z^z) and [z^z). The kernel K^^''^^ [z^z] was denoted by 

Pr{x)6{y) + pc{z)/2 in Refs. [HI El] which is the sum of the distribution of the right 
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handed modes and the half of the distribution of the complex eigenvalues. Then the 
kernel {x, x) is equal to the chirality distribution over the real eigenvalues 

p^{x), see Refs. [22l HH El] . 

Again the Pfaffian determinant as well as the identification with two flavour 
partition functions of carry over to the microscopic limit and, thus, to eigenvalue 
correlations of the Wilson-Dirac operator in the e-regime. The microscopic limit of 
the three kernels k[~^''^\ K^^''^^ and K^'^''^^ are trivial corollaries of Eq. fl3.59l) . The 
derivation of this limit for the other three kernels will be made elsewhere [63] . Also the 
discussion of the results will not be done here. 

Again we can look what our results mean inn the notation of Refs. pOf [211 [22] . In 
the microscopic limit the following kernels are proportional to the two flavour partition 
functions of the chiral Lagrangian 

Ak[~^''^\zi, Z2) oc (^1 - X2) Im Z^Jrne -xi,mQ- X2; A7, A7; as, ag/r = 0)5{yi)5{y2), 

£1— >o ' 

£2^0 

(4.33) 

AK'^ '{zi,Z2)(xlm d{y2), (4.34) 

e-s>0 X2 — Z\ 

^6"^'°^ ''(^15^2) OC (Si - 22)^2/0(^6 - ^1,^716 - ^2; A7,A7; 08,06/7 = 0). (4.35) 

The other kernels are only integral transforms of these three partition functions. As 
for we can create the case 00/7 7^ by multiplying the expression fl4.18p with the 
partition function of A^f fermionic flavours and integrating over Gaussian distributions 
of mg and A7. At the end we divide the resulting expression by the partition function 
with Af fermionic flavour and with 05/7 7^ 0. 



5. Conclusions 



We derived the orthogonal and skew-orthogonal polynomials corresponding to the 
Hermitian as well as non-Hermitian Wilson random matrix ensemble. The orthogonal 
polynomials are the Hermite polynomials from order to — 1 in both cases. They 
result from the i/-dimensional GUE describing the broadening of the v generic real modes 
which are at zero lattice spacing the zero modes. Such a GUE was already discovered 
in the chirality distribution over the real eigenvalues ESI 121] as well as in the level 
density of the Hermitian Wilson random matrix ensemble and, thus, the Wilson Dirac 
operator [201 ED (32]- Surprisingly this GUE is already the universal result and is a 
dominant part in the eigenvalue correlations at small lattice spacing since it forms the 
Dirac delta functions at zero with weight v in the continuum limit, see Refs. [211 [22] • 

The remaining spectrum is described by skew-orthogonal polynomials starting from 
order v. They describe the remaining spectrum apart from the v generic real modes. 
In a unifying way we constructed these polynomials and derived recursion relations 
which enable us to obtain the odd polynomials by simply acting with a derivative 
operator on the even ones, cf. Eq. (13.311) . This derivative operator can be identified by a 
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creation operator of a harmonic oscillator. Moreover we derived a Christoffel Darboux- 
like formula f l3.37p which is equivalent to the partition function of two fermionic flavours, 
see Eq. (13.531) . The even skew-orthogonal polynomials are equal to one flavour partition 
functions, see Eq. fl3.43l) . With help of this knowledge we were able to derive the 
Rodrigues formula (13.451) interpolating between the one of the Laguerre polynomials 
and the one of the Hermite polynomials. 

As an application we considered the unquenched fc-point correlation functions of 
the Hermitian and non-Hermitian Wilson random matrix ensemble. We derived a 
Pfafiian factorization in both cases. The one of the Hermitian matrix was already 
know before [23] but we traced the entries back to the two flavour partition functions, 
see subsection 14. 1^ which is numerically advantageous. The Pfaffian of the non- 
Hermitian random matrix is a completely new result. We identified its kernels as 
two flavour partition functions, too, see subsection 14.21 These partition functions 
can be readily interpreted as correlations of the complex conjugated pairs, the real 
eigenvalues corresponding to the right handed modes and the average chirality over the 
real eigenvalues. 

Although the random matrix is non-Hermitian we did not need a Hermitization 
as it was introduced in Ref . [65] . We circumvented this approach by splitting the kernels 
with bosonic flavors into two kinds of terms. One kind corresponds to the chirality over 
the real eigenvalues which exhibits no singularities in the bosonic determinants. The 
other term are integral transforms of partition functions with fermionic flavours instead 
of bosonic ones. Hence there are no problems of integrability anymore. Especially 
we have not to double the number of the bosonic dimensions in the superspace when 
applying the supersymmetry method. 

The Pfaffian factorization as well as the identification with two flavour partition 
functions carry over to the microscopic limit and, thus, to the spectral properties of 
the Wilson-Dirac operator in the e-regime [271 EH EH [30]. Hence the results shown 
in Sec. S] are a good starting point for an analytical study of the Hermitian and 
non-Hermitian Wilson-Dirac operator. In particular the calculation of the individual 
eigenvalue distributions will benefit of the structure since a representation as Fredholm 
Pfaffians are possible, see Ref. [61]. Fredholm determinants and Pfaffians are compact 
expressions simplifying the perturbative expansion of the gap probability in the /c-point 
correlation function to obtain the individual eigenvalue distributions. 

Moreover, the skew-orthogonal polynomials and the Christoffel Darboux-like 
formula also appearing as kernels of the Pfaffian determinants reduce to a quickly 
converging sum, see Eq. (I3.58p . and two-fold integrals over phases, see Eqs. (I3.59P and 
(I3.60p . respectively. 

The Pfaffian of the A;-point correlation function will persist in the continuum limit. 
It is in agreement with Ref. [62j where a non-trivial Pfaffian determinant was derived 
for (3 = 2 random matrix ensembles. A similar but not completely equivalent structure 
was derived in Refs. [661 EZ] for /3 = 2 ensembles, too. Hence Pfaffians seem to be more 
universal than the determinantal structures in the eigenvalues statistics of RMT. 
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The Pfaffian determinants we found reflect the breaking of the generic pairing of 
eigenvalues in the continuum limit to no reflection symmetry at all in the Hermitian 
case and the reflection symmetry at the real axis in the non- Hermitian one. For example 
at a = and fixed matrix we can say if A is an eigenvalue of then —A is also one. 
This is not anymore true at finite a. 
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Appendix A. De Bruijn-like integration theorems 



We generalize the de Bruijn-like integration theorem [6^ to an integrand which 
is a product of one determinant and one Pfaffian, see Appendix A. 1, and of two 



determinants, see [Appendix A. 2 



Appendix A.l. With a Pfaffian integrand 



Let Ai, A2 and A3 be three positive integers fulfilling the condition 2 A3, A2 > Ai > 0. 
We consider the following integral 
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Pf 
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{-Dt{zc)} {Etc} 
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d[z]. 



(A.l) 



The matrix B is an arbitrary constant matrix whereas E is an antisymmetric constant 
matrix. The matrix valued functions A, C and D are sufficiently integrable and C is 
antisymmetric in its entries. 

After an expansion of the determinant in Eq. (lA.ip in the entries Ac{zb) we can 
integrate over the variables z [57], i.e. 

Ii = TTr—i:TT7 ^ signa;det[5fe^(c)] (A.2) 
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The remaining determinant can be combined with the Pfaffian by the sum. Thereby we 
use the identity 



det[5,^(,)] = ^-l)(N.-m)iN,-N.-l)/2 

Ni+l<c<N2 



X Pf 







{Bcu{b)} 



l<b<N2-Ni 
Ni+l<c<N2 



Ni+l<b<N2 
l<c<N2-Ni 





This yields the result 



= (-1 


■^iN2-Ni)iNi+N2-l)/2 j<^^^ 








j Ab{z^)A,{z2)C{zi,Z2)d[z\ 


1 Ak{z)D,iz)d[z] 


Bcb 


X Pf 


- 1 A^^,){z)Dh{z)d[z] 


Ebc 







—Bbc 









(A.3) 



(A.4) 



The number of the first set of coloumns and rows is N2, the one of the second set is 
2A''3 — A^i and the one of the third part N2 — Ni. Hence, we take the Pfaffian of a 
2{N2 + N3- Ni) X 2{N2 + N3- Ni) antisymmetric matrix. 



Appendix A. 2. With a determinantal integrand 

Now we study the integral with a determinant instead of a Pfaffian, cf. Eq. (lA.ip . i.e. 

{Ac{ZbR)} 



det 



{B,{ZbL)} 



l<b<Nii 
l<c<Ni 

l<b<NL 
l<c<Ni 



{Cbc} 



l<b<Ni-NR-NL 
l<c<Ni 



det 



{D{ZbR,ZcL)} {E^iZbR}} 

l<b<NR l<b<NR 

1<c<Nl 1<c<N2~Nl 

{Fb{z,L)} {Hbc} 

l<b<N2-NR l<b<N2-NR 
1<c<Nl 1<c<N2-Nl 



d[z] 



(A.5) 



The matrices C and H are arbitrary constant matrices and the matrix valued functions 
A, B, D, E and F are chosen such that the integrals exist. The positive integers Ni, 
N2, Nn and Nl have the relations Ni> Nr + Nl and N2> Nn,NL. Without loss of 
generality we can assume A'"^ > Nr > 0. 

In the first step we split both matrices F and H into two blocks, i.e. 

{Fbiizci)} 



[FbiZcL)] 



l<b<N2-NR 
1<c<Nl 



{Fb2{ZcL)} 



l<b<NL-NR 
1<c<Nl 

l<b<N2-NL 
1<c<Nl 



bc\ 



l<b<N2-NR 
1<c<N2~Nl 



{Hbci} 



l<b<NL-NR 
1<c<N2-Nl 



{Hbc2} 



l<b,c<N2-NL 



(A.6) 



(A.7) 
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where we assume that H2 is invertible. Later on we will relax this restriction since 
I2 is a polynomial in the constant matrices C and H. We pull H2 out of the second 
determinant and have 

{AciZbR)} 



I2 = det H2 / det 



l<b<Nji 
l<c<Ni 

l<b<NL 
l<c<Ni 



l<b<Ni~NR~NL 
l<c<Ni 



(A.8) 



X det 



DizbR, z^l) - ^ Ei{zm) [H^ ^) ^- FjiizcL) ^^^^^^ 



FbliZcL) - 



l<i,j<N2-NL ) 



<b<NL-Nii 
1<c<Nl 



d[z]. 



After an expansion in both determinants we obtain 

Nr 

' ^^ic{b+Na){ZL) 



Nr „ 

l2 = NL\detH2 signwJJ / 



(A.9) 



X 



D{ZR, zl) - ^^(^«) (^2"'),,- Fi2{zl) 

l<i,j<N2-NL 

Nl-Nr 



d[zR, zl] 



X 



F,,{zl)- Y Hui{H^')^^F,2{zl) 

l<i,j<N2-NL 



d[zL] 



Y\_ / B^ib+2NR,){ZL) 
b=l 

Ni-Nr-Nl 

X n Cbu;{b+NR+NL) ■ 

6=1 

Notice that the sum over the permution of the second determinant can be absorbed into 
the first one which gives A^^^! . 

To shorten the notation we define the folowing matrices which are integrals over 
one or two variables 



Obc = j {Ab{zR)Bc{zL) - Ac{zR)Bb{zL)) D{zr, ZL)d[zR, zl] , 



be 



Ab{zR)Ec{zii)d[zR] 



Qbca= J FbaizL)B^{zL)d[zL] , a G {1, 2} 
Then the integral ( 1A.9P reads 
I2 = Nl\2~^'' det H2 Y siS^'^ 



(A.IO) 
(A.ll) 
(A.12) 

(A.13) 



X 



Nr 

n 

b=l 



Ouj{b)uj{b+NR) — 'Y^ {Pui{b)iQ jui{b+NR)2 — Pui{b+NR)i 
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Nl-Nr 








Ni-Nr- 


Nl 


X 


n 


Qbuj{b+2NR)i — 


Hbii {H2 ^)^■ 


Qjijj{b+2NR)2 


n 


Cbuiib+NR+NL) ■ 




b=l 


l<i,j<Af2- 






b=i 




This 


sum can be represented as a 


Pfaffian, i.e. 








h = 


I^_^^Ni{Ni-1)/2+Nr{Nr+1)/2j^^ 


Nr\ det H2 Pf 


Rbc 


Ccb 




(A.14) 










-Cbc 








with 

[-Rbc] 



Obc Qcbl 
-Qbcl 



(A.15) 



+ 



E 



Qib2 




Pbi 
-Hbil 



{H,-')„ 



Qjc2 





-Hcji 



Pushing the determinant of H2 into the Pfaffian we have the final result 



X Pf 



X Pf 



Ni{Ni~l) 
Obc 


/2+Nr(i 
Qcbl 


Vj? + 1)/ 

Ccb 


'2+{N2- 
Qcb2 


-Nl){N2-1 
Pbc 


—Qbcl 











—Hbci 


—Cbc 














—Qbc2 











—Hbc2 


—Pcb 

Ni{Ni-l) 

Obc 


Hcb 

/2+Ni 
Qcb 


1 




Vh + 1)/ 
Pbc 


Hcb2 

I2+{N2' 

Ccb 




-Nl){N2-1 


—Qbc 





—Hbc 





—Pcb 


Hcb 








—Cbc 












(A. 16) 



The dimensions of rows and columns are from top to bottom and left to right 
{Ni,N2 - Nr,N2 - Nl,Ni - Nr- Nl). In Eq. ([06|) we drop the invertibility of 
the matrix H2 because I2 is a polynomial of this matrix. 

Appendix B. Derivation of the coefficients in the recursion relation 



In Appendix B.l we show that the recursion relations of the polynomials g^^^ take the 



form f l3.3ip and f l3.32p . The coefficients e[^'' and ej^"* are derived in Appendix B.2 
Appendix B.l. The general form 



In the ffist step we take the scalar product fl2.25p of Eq. f l3.3Up with p^i^\ < k < u — 1. 
We find 



/q(±) _ ^ /n(±)„(±)|J±)\ - /nW|r,W'\ I /^r T/il , (±)| (±)^ 
Plk - %+l\Pk /.(±) - - — \%+l\Pk /„(±) \1u+l\Pk /g(±) 



hi 



9l 



hi 



' 9{ 



2hk 



0, (B.l) 
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Thereby we used the fact that gf,^] and p]^' , < k < u — 1, are orthogonal to each 
other. 

The monic normahzation of the polynomials enforce the condition 



a 



(±) 



n 



(B.2) 

For the other conditions we take the anti-symmetric product of Eq. f l3.30p with q'j^ii- 
Let / = 2i and /' = 2i' we have 



Hu+2i\%+2ii J (±) ~ ydu+2i\^ %+2i' 



92 

_ (i) (±) p> _ 

where the integrated Kronecker delta is 



(±) (i) p 



1, /eNo, 

0, else. 



Combining Eqs. f lB.2p and flB.3p we find 



'^2i,2j'+l ~ '^2i',2j+l ~ a2 



(B.3) 



(B.4) 



(B.5) 



With I = 2i + 1 and /' = 2i' + 1 we get another relation 



n(±)o(±) |o(±) ^ _ _ (J^) |n(±)oW 

'du+2i+l\%+2i'+l (±) ~ iyj^+2i+ll-^ yi/+2j' 
''92 ^ 



+1 / (±) 

^9^2 



(±) (i) p. 



2j'+l,2j®«'-«+l' 



(B.6) 



With this we conclude 



„(±) 



'^2i+l,2j' 



n o. 



(±) 



2 ^(±) 



a 



0, 



i-i 



i + 1, 



(B.7) 



0, else. 

The constants e-^^ cannot be specified by Eq. (IB. 61) . 

The last relation which we get by the skew-orthogonality of the polynomials is the 
one for the choice / = 2i and /' = 2i' + 1, i.e. 

■ W in(±)^W 



^ 'iu+2i\%+2i'+l I (±) 
92 



)(±)^(±) u(±) 
The identity yields 



Oi' Oi2i,2i'^i-i' 



Hu+2i\-^ Hu+2i'+l 
i "2ii'+l,2i+l^*'-«' 



„(±) 



2i,2j' '^2j+l,2i' + l 



(B.8) 
(B.9) 

Again the constants ej^^ have to be determined. 

Collecting the intermediate results (IB.ip . (IB.Sp . ( 1B.7P and ( 1B.9P the expan- 
sion fl3.30p reduces to the results f l3.3ip and (13.32 p . The derivation of the constants 
e] and remams. 
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Considering the anti-symmetric product of Eq. f l3.3ip with pf'^21+1 ^"^^ 
\ 
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u+21+1 1 (±) 

92 



yu+2i+i I (±) 



92 



(±) 

!^+2'+l / {±) 



(±) 



f'T T f'l (±) 

o Piy+2l+l 



(±) iJ±) 
u+21+2) j±) 

92 



I 

o(±)^±) 
^(±)^±) 



(B.IO) 



A similar calculation can be done with the scalar product of Eq. fl3.3ip with pl^2i 



9i 



(z/ + 2l)K+2l-l 



i/ + 2«l/^i/ + 2«-l)g(±' 



yUr =F yUl 



l^v+2l^l 

r^u+2i^i 
r^v+2i^i ■ 



(B.ll) 



In both calculations we employed the definition (13. 5p . the orthogonality of the 
polynomials pj^^ and Eqs. ( KWf . ^72E\f . ^72E\f and (1^:^ . The combination of both 
results yields the recursion relation 



%+2l \Pu+2l+2 



f2 



n 



-(/ir =F/ii) + 



9'!'+2i-2bl+2i 



(B.12) 



The starting point of this recursion is Z = 0. Due to the definition (13. 5 p we know that 
pi 



„(±) 



""■^^ and ql^i = p^^i- We conclude 



2n 



2n 



(±) 



Pu+l 



a o. 



n 



Hence, we can solve the recursion and find 



%+2l\Pu+2l+2 j 



(l + l)a'o[^^ 



(yWr =F/ii), 



{21 + 1) 



(B.13) 

(B.14) 
(B.15) 
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In a similar way we derive an identity for the constants e;^"*. We take the scalar product 
of Eq. (13.321) with p[^2i obtain 



9i 



(±) 



h 



u+21 



(±) \J±)'\ , /^r , (±) I (±) V 

{Qu+2l+l\P,y+2l) g[±) 



iv+2l+l 



+ 



n 



1 



(±) 



(±) lj±) 

/^r =F /^l / (±) I (±) \ \ 

/ (i) I (±) ^ 



u+21 



{iy + 2l){q]f^>2i+i\Pu+2i-i)g[^) 



n 



n 



(z/ + 2/ + l)a2 2 



/ (±) I (±) ^ 



+ 2/ + l)a2 <^^^^^SW2bSW3(±) 



9i 



v+21+1 



h 



v+21-1 



(/ + T^l) , (±) I (±) . 



■2«+l 



I h^+2i+i 
{I + l)(/ir T/Ui) 



/ (i) I (i) ^ 

\1u+2l+l\Pu+2l-l) a^±) 



9i 



i'u+21-1 



„(±) |J±) 
yi/+2« l/^i/+2«+2 / (±) 



n 

7^ 



'I 

I (±) I (±) ^ 
\%+2l+i\Pu+2l+lI gf) 



92 



/ (i) I (±) ^ 
\1u+2l+l\Pu+2l-l) r.^±) 



9i 



+2«+l 



+2«-l 



n 



(B.16) 



For a further simplification we need more information, i.e. we have to perform the 
integral {qlX\i+i\pf+2i-i) g^^^ ^^t all / e Nq. 



Appendix C. Derivation of the Christoffel Darboux-like formula 



Let Zi and Z2 be restricted to the real axis, i.e. 2:1/2 = Xi/2- The action of the sum of 
the two differential operators D^^'^ with respect to xi and X2 on S^'^^ is 



d d 
+ 



n 



dxi dx2 



{Xl + X2) + (/ir ± yUl) ) T!~^\{xi,X2) 



(C.l) 



1=0 



(±) i^^^det 



n 



+ 



(±) 

n—l 



no. 



(±) 



a^o 



(±) 
i-i 



n , 
+ ^det 



det 



(iS2l+2i.^2) q\S2li.^2) 

q^^+2i^2i.^i) ^SU^^i) 



det 



(C.2) 
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Let X = (xi + X2)/2 and Ax = {xi — X2)/2. Then we rewrite the differential equation 
to 



dX 



exp 



n 



- ^X^ + [fi, ±fii)x] Eit\ (X + Ax, X - Ax) 



(C.3) 



(±) 

n-l 



exp 



n 



--X' + {fi, ±/ii)X 



det 



qit\n-2iX - Ax) ql^UX-Ax) 



In the next step we itegrate this equation from X to cxd and take into account that the 
upper boundary vanishes due to the Gaussian and yields Eq. f l3.37p for real entries. The 
restriction to real Zi and Z2 can be relaxed since the integrand is absolutly integrable 



Appendix D. Derivation of Eq. (I3.44p 

We consider Eq. (13. 43 p . The characteristic polynomial in zf''^''^^ can be raised into the 
exponent by a Gaussian integral over a complex vector of Grassmann (anti-commuting) 
variables, 



6 

6 



Ir 



61 



Ir ) 



Sir ; ?ll ■■ 



(D.l) 



The integration is defined by 
^d^, = [ 



and / di^ = / d^* = 



(D.2) 



Moreover we employ the conjugation of the second kind, i.e. 

(O* = -e. and (6e,)* = Ce;. (D.3) 

Good introductions in the standard techniques of supersymmetry can be found in 
Refs. 



We find 

oc / exp 



n 



Id- 



(tr + tr B^) - n tr WW'^ + tr A{ii^ + i\ 



(D.4) 



X exp 



tr i?(/ii + 6e/) + tr Wixil - tr W'^i,^, + z{i\i, ± e/^i )! 



oc / exp 



|^(tr(/ir + 6e/)' + tr(/xi + 6e/)') - \ 



n 



X exp 
oc J exp 

X exp 



1 



n 



n 



Z ) ilir 



n 
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With help of the superbosonization formula [52| 153] we express the integral over by a 
two-fold integral over two phases (I3.44p . 

Appendix E. Derivation of Eq. ( I3.53h 

Equations f l2.27p and fl2.45p with the parameters k = = h = and iVf = 2 read 



{-Zi, -Z2) OC 



Zi — Z2 



(E.l) 



X Pf 









hiSij 




















'?t.+2i+l(^l) 


^1^21(^2) 
ql%i+i{z2) 




















Hv+2j+l 













. -Pf\z2) 


















combined with the derived knowledge in subsection 13.11 This Pfaffian can be expanded 
in the normalization constants hj. Then the polynomials p\'^\ < I < drop out. 
Furthermore, we make use of identity (14. 3 p and get 

^ (E.2) 



Z. 
X Pf 



2 i-Zi,-Z2)(X 



Zl - Z2 



Sr^J_i Tu+2M^) Tu+2j+li.Zl) 



-1 

1 



qu+2i (^1) qlZ'2i{^2) 



(±) 

qv+2i+l qv+2i+l 



—(+S (^!'+2i(^l)^^+2j+l(^2) - 5'!-t2j+l(^l)^^+2j(^2) 

Zl — Z2 . o- ^ 

This result is proportional to the sum cf. fl3.35p . 



(E.3) 



Appendix F. Simplification of the kernels 



In Appendix F.l and Appendix F.2 we simplify the kernels of D5. Derivations of the 



kernels of are given in [Appendix F.3[ [Appendix F.4| and [Appendix F.5 



Appendix F.l. The kernel k[ '""^ 

With help of Eqs. ^^MM) it can be readily shown that 





1 


u—l _ n—1 

n ,n 

i=o ^ i=o 


1 






2n+v 


(5^) 








(2n + z/)! 




I 










r 'A 


{Xi , Xj ) 






Xi) 


'A 


\xi,X2) 






'\^^) 


X Pf 




'\xi,Xj) 
'\x2,Xj) 


A- 




^(3^2, 


Xi) 


A- 


\xi,X2) 




xP{ 
X2 


-'Ac 


\X2) 








-4- 


-'si- 


\xi) 


— X2 


''9['\X2) 






(F.l) 
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The indices i and j run from 1 to 2r;, + in the first row and column and 1 to z/ in the 
last ones. 

In the next step we extend the Vandermonde determinant by two Dirac delta 
functions such that the integration is over 2n + z/ + 2 variables, 



< ' n^n 



(2n + u + 2)\ J-A/i. 11 (-) 



d[x\ 



(F.2) 



X det 



X 



i-i 



6{xj — xi) 
6{xj — X2) 



Pf 



■ 'A- 




xi g\ '[xi) 


—Xj 








The two Dirac delta functions can be expressed by the imaginary parts of the Cauchy 
transforms in two variables. Using the identity 



det 



X, 



1 



Xj — Xi — ISi 



^ -|^^(2n+i/+2)(2n+i/+l)/2 {^l + ISl - X2 - 162) ^2n+u+2{x) ^ 



Xj — X2 — '1^2 

we find the expression 



2n+v+2 

n {xj - xi - iei){xj - X2 - ie2) 
i=i 



Kr"\x,,x2) 



u—1 ^ n— 1 



1 X2 — Xi 



(2n + z/ + 2)! J-_l/,^. -_- ^ 



Im 



j=0 j=0 j £2-i>0lR2n+i'+2 



(F.4) 



A 



X 



2n+u+2[X) 



2n+u+2 

Yl (Xj — Xi — iei){Xj — X2 — 162 



-Pf 
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xi g\ '{Xi) 








This result is the partition function of with two bosonic fiavours, see Eq. (14. 8p . 



Appendix F.2. The kernel K2 '"^ 

Again we start from an identity between the kernel and an integral weighted by the 
joint probability density ps, see Eq. (12. 6p . i.e. 

u—l ^ n— 1 



Xi,X2j 



(2n + 



1 n-n-- 



d\x\/^2n+u{.X,Xi] 



(F.5) 





r 9^ 






9i- 


''(a;i,X2) 


{-) 
x^ ^1 




X Pf 




\X2^ 







i-i (-) 
A 9i 








-'9!- 




X2 


"^^71"^ (^2) 








Notice that we integrate this time over 2n + u — 1 variables. Hence the range of the 
indices i and j is from 1 to 2r;, + z/ — 1 in the first row and column and from 1 to z/ in 
the last ones. 
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The integral is extended to 2n + v variables by introducing a Dirac delta function, 

^_lY2n+i'-l){2n+u-2) /2 



K[ '''\xi,X2 



{2n + 1^ + 2)1 



U—l 11— 1 

n^n 







1 

FT 



d\x\ 



2n+u 

n (^i - ; 



X2 — Xi 



(F.6) 



X det 



Pf 





x\ ^g[ \xi) 








5{xj - X2) 

We employ again the Cauchy integral as a representation of the Dirac delta function 
and an equation similar to Eq. fIF.Sp . This yields the result f l4.9p which is the partition 
function of with one fermionic flavour and one bosonic one. 

Appendix F.3. The kernel /Tg'*''"'* 

Also for the kernels have a representation as an integral over the eigenvalues 
weighted by the joint probability density ( !2.14p . 

/I \n(n+l)/2+!/(!/+l)/2 '^^l " ^^"^ 



K. 



n\{n + z/)! 



j=0 J j=0 . 



(F.7) 









92 [^i 1^1) 


X det 






9''^\z2,zi) 








-x\-'g\^\x,)5{y,) 



The index j in the first column takes the values 1 to n + z/ while i goes from 1 to n in 
the first row and from 1 to in the last one. 

We expand the determinant in the row with the variable Z2 and have 

K^^'''\z,,Z2)=g^^\z2,z,) (F.8) 



+ 



(^_iyi(n+l)/2+{v+2){u+l)/2 

n\{n + V — 1)\ 

9^'i^\4') 



u—l n— 1 



d[z, Z\9^2^^ {Z2, z)A2n+y{z, z) 



X det 



92 V'i ?^iJ 



The integration of the second term is extended by a Dirac delta function. However this 
distribution can only be symmetrized with respect to the l!^^ integration in contrast to 



the calculation in Appendix F.l and Appendix F.2 We add and subtract a Dirac delta 
function for the integration over 'zS'^\ Collecting these steps we find 

qt\z2.zi) (F.9) 



^^"H^l,^2 



+ 



(ra - l)!(n + z/)! J-i/i — 



d[z, z]g^^\z2, z)A2n+u{z, z) 



X det 



9'2'\zhz)'') 
9i'\t\f) 
i^fr'9^'\xf)5{y^P) 



<l)^ 
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d[z, z\ — - — —dei[z - zl2n+v) 



2n + i/ + l 





-5^^)W-zl)' 







det 



y2 i^j ? 



— X 



r(l)^ 



is needful to construct the chirality 



with = 6{x)5{y). 

The minus sign in front of S^'^^zf^'' 
distribution over the real eigenvalues. Thereby we need the following relation 
between the real eigenvalues of +?7?.l2n+i/, A^-^^+m, and the eigenvalues of +77175 

m^M, (F.IO) 



dm 



\(5) \(W) , „ 
K =K +m=0 



(W) 



(F.ll) 



where ipi is the eigenvector to the eigenvalue A ■* of • The right hand side of 
Eq. (IF.lOp is the chirality of the corresponding eigenvector. Since the eigenvectors of 
the complex eigenvalues have vanishing chirality Eq. fIF.lOp is only applicable to the 
real modes of . The following short calculation will show the connection between 
the chiral distribution over the real eigenvalues and the third term in Eq. (IF.Qp . 
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Hereby we have to understand the whole calculation, in particular the limit of the 
imaginary increment le, in a weak sense. The complex conjugated pairs of the integration 
variables in the third term of Eq. f lF.Qp do not contribute. We recognize this by 
expanding the second determinant in the two point weights and the first determiant 
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If 2i 



z the integrand is anti-symmetric under the complex conjugation of z, 



i.e. z o z*. Due to _ = 5(2)(;^* _ 

the determinant is symmetric under 
z z* while (yf2'''''(-2;, -2*) ~ gdz) = —gdz*), see Eq. f l2.17p . Thus the integral over 
the imaginary part of z vanishes. The same discussion can be made for all complex 
conjugated pairs. 

Expanding the determinant in the Dirac delta function and using the 
calculation (1F.12P we find 
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The second term is an integral transform of the partition of with two fermionic 
flavours and the last term is an integral over the partition function with one bosonic 
and one fermionic flavour. Notice that the integral over z' does not commute with the 
limit £ — 7- because of the singularity at Xi. This singularity cancels with a term after 
we take the limit. Hence, the expression (1F.14P is equal to the result fl4.30p . 



Appendix F.4- The kernel K\ 
The starting point for this kernel is the identity 
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The index j runs from 1 to + z/ — 1 and the index i takes the values 1 to n + 1 in the 
upper row and from 1 to z/ in the lower one. We introduce two Dirac delta functions 
and, thus, extend the integral by two additional i^'^ variables. 
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We extend the determinant by Dirac delta functions of z} similar to the calculation in 
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In the final step we expand the last three terms in (72^''(z*, i^^'') and the Dirac delta 
functions can be rewritten as limits of Cauchy transforms, see Eq. flF.12p . 
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Only the first term is new in comparison to the kernel i^^^"*"'"^ and it is the partition 
function of with two bosonic flavours which agrees with the result fl4.3ip . 



Appendix F. 5. The kernel 
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where the indices of the determinant are j G {1, . . . , n + z/ — 1}, i G {1, . . . , n} in the 
first row and i G {1, . . . , i/} in the last one. The extension with a Dirac delta function 
yields 
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We proceed in the same way as in [Appendix F.3| by extending the first determinant by 



^1) 



Z2) and expanding the resulting correction in g^\z2,zf^). Then we find 
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The first term is the partition function of Dy/ with one bosonic and one fermionic flavour 
and the second term is the kernel K^'^'"\ Therefore Eq. flF.2ip is the result fl4.32p . 
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